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Abstract 

Precise energies of rovibrational states of the exotic hydrogen-like molecule {dt^)Xee are of 
importance for dt^ resonant formation, which is a key process in the muon-catalyzed fusion cycle. 
The effect of the internal structure and motion of the dtfi quasi-nucleus on energy levels is studied 
using the three-body description of the {dtfi)Xee molecule based on the hierarchy of scales and 
corresponding energies of its constituent subsystems. For a number of rovibrational states of 
{dtn)dee and {dtfi)tee, the shifts and splittings of energy levels are calculated in the second order 
of the perturbation theory. 

PACS numbers: 36.10.-k, 36.10.Dr, 33.20. Wr 
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I. INTRODUCTION 



It is known that one stopped muon in a deuterium-tritium mixture yields more than 100 
nuclear fusion reactions. The process of muon-catalyzed fusion has been intensively studied 
and a detailed description exists in the literature, e. g., in review articles 

mm 

One of 

the key processes in the muon-catalyzed fusion cycle is the formation of the hydrogen-like 
exotic molecule {dtfi)Xee (for the sake of generality X stands for either isotope d, t, or p), 
in which a dtfi mesic molecule substitutes for one of the nuclei in the hydrogen molecule. 
It is widely accepted that the resonance mechanism proposed by Vesman ^ is responsible 
for the high rate of the dtfi formation. Due to this mechanism, a dtfi mesic molecule in a 
loosely bound excited state to be produced by low-energy collisions of tfi mesic atoms and 
DX molecules in a resonance process 

tn + DX — > {dtn)Xee 

followed by dtfi transition to the tfi ground state. The rate of the resonance process is 
sensitive to the precise resonance position and an accuracy better a 1 meV is necessary to 
obtain reasonable theoretical estimates of the formation rate 0,0,0, ^. 

Resonance formation can take place if the energy released in dtfi binding is transferred 
to the rovibrational excitation of the exotic molecule {dt^)Xee. This is actually the case as 
dt^ has a loosely bound excited state with an angular momentum A = 1 and binding energy 
which is comparable to vibrational quantum of the {dtfi)Xee molecule. In a non-relativistic 
approximation, different calculations determine with a good accuracy the binding energy 
of the isolated dtfi mesic molecule . To obtain the precise value of the binding 

energy one has to correct the non-relativistic energy for relativistic effects, hyperfine effects, 
finite nuclear size, vacuum polarization, and others. The resonance position is determined, 
besides the binding energy of isolated dtfi, by the energy of the rovibrational excitation of 
the hydrogen-like molecule {dt^)Xee with one nucleus being the particle X and the other the 
excited dtfi mesic molecule. As the "size" of the excited dtfi mesic molecule with A = 1 is of 
the order of 0.05 a.u. 0, which is much smaller than the internuclear distance in the whole 
molecule, the rovibrational spectrum of (dtii)Xee can be calculated to a good approximation 

nn 

by treating dtfi as a point-like charged particle p, 17( . Nevertheless, to reach an accuracy of 
the order of a tenth of a meV one should take into account the energy shift which arises due 
to the internal structure and motion of a dtfi mesic molecule. 



The effect of the dtfi finite size was previously investigated in a simple approach 
y, ^1 where the energy shifts for the {dtfi)dee were obtained by multiplying by 1.45 the 
shift calculated for the atom-like system {dtfi)e in the second order perturbation theory 
(PT). Within the framework of this simple approach it is not possible to take account of 
the molecular structure; in particular, the calculated energy shift is independent of the 
rovibrational quantum numbers. The effect of the molecular structure, i. e., the dependence 
on angular momentum, was explicitly demonstrated in the elaborate six-body calculation ^ 
of the {dtii)dee energy shifts in the first order of the perturbation theory. Note, however, that 
the first- and second-order PT contributions to the energy shift are comparable. Recently, 
resonance positions in the low-energy t/i + D2 scattering have been obtained in the elaborate 
three-body calculation 

HQ. 

Only few resonance states with the zero total angular 
momentum have been considered in this paper. 

The main aim of the present paper is to calculate the energy shifts which arise due to 
the internal structure and motion of the dt^ mesic molecule embedded in the hydrogen-like 
{dtii)Xee molecule. The calculation is reduced to solution of a three-body problem for heavy 
particles t/i, d and X. This approach is based on the hierarchy of scales and corresponding 
energies of constituent subsystems of the {dtfi)Xee thus reliably taking into account the 
specific features of this molecule. As a result, the energy shifts are obtained for a number 
of vibrational and rotational states of {dtfi)dee and {dtfi)tee in the second-order PT. 



II. METHOD 



The structure of the exotic molecule {dtfi)Xee is characterized by a hierarchy of scales 
and corresponding energies of its constituent subsystems. In this respect, a tfi mesic atom 
is small in comparison with its mean separation from a deuteron in the loosely bound dtn 
mesic molecule, which allows to be treated as a point-like neutral particle interacting 
with a deuteron by the short-range effective potential. There is also interaction of tfi with 
the second nucleus X; however, this might be neglected due to large separation between 
these particles. In turn, the size of a dtfi mesic molecule is small in comparison with the 
amplitude of vibrations in {dtfi)Xee] therefore, it moves as a point-like quasi-nucleus near 
the equilibrium position. For this reason, the effect of the dtfi structure is considered within 
the framework of the perturbation theory. 
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Furthermore, two electrons in the hydrogen-hke molecule {dt^)Xee move much faster 
than the heavy particles d, X, and tfi, which makes it possible to use the familiar Born- 
Oppenheimer (BO) approximation, i. e., to solve electronic problem with the fixed charged 
particles d and X thus obtaining the BO energy which plays a role of the effective poten- 
tial between d and X. The electronic excitations, which require a considerable amount of 
energy are not taken into account for the low-energy processes under consideration. 

As a result, the description of {dtfi)Xee is reduced to solution of a three-body problem 
for three particles tfi, d and X. The interaction between charged d and X is described by 
the well-known BO potential for the hydrogen molecule. In accord with the treatment of 
the t/i mesic atom as a point-like neutral particle, the present calculation does not explicitly 
use the tfi + d effective potential, rather the result is expressed via the low-energy tjj, + d 
scattering phase shifts and characteristics of the mesic molecule in the loosely bound 
excited state. 

The {dtfi)Xee states are either true bound states or narrow resonances if their energy is 
below or above the tfi + DX threshold. As the energy shifts are mainly determined by the 
coupling with closed channels, in the present calculation both resonances and bound states 
are treated on an equal footing thus neglecting a small contribution to the energy shifts 
which comes from the coupling with the open t/i -|- DX channel. 



A. Three-body description 

Under the above approximations, the Schrodinger equation for the hydrogen-like molecule 
[dtii)Xee reads 



1 A 1 
A,. 

2/ii 2/i2 



A, - —A, + Vi{r) + V2{\p-(3v\) + V{\p + av\) - E 



^ = (1) 



where the Jacobi coordinates r and p are the vectors from d to the point-like mesic atom 

t/i and from the second nucleus X to the dt^ center of mass, respectively. The reduced 

mim2 (mi + 1122)1713 mi 
masses and parameters a and p are /xi = , /i2 = , a — 



mi + m2 mi + m2 + m.3 mi -|- m2 

and (3 = , where mi,m2 and m^ are the masses of t/i, d, and X, respectively. The 

mi + m2 

atomic units are used throughout the paper unless other is specified. In Eq. ((T}, l^(|p + Q;r|) 
denotes the well-known BO potential describing the interaction between charged d and X 
while the short-range potentials Vi(r) and V2{\p — /3r|) describe the interaction of a tfi 
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mesic atom with a deuteron and X, respectively. In the following, due to large internuclear 
separation (p ^ r) in {dtfi)Xee, the short-range interaction V2(|p — /3r|) of the tfi mesic 
atom with the second nucleus X is negligible and will be omitted. 

A natural zeroth-order approximation for the calculation of the {dtfi)Xee energy levels 
is to treat the dtfi mesic molecule as a point quasi-nucleus with the dtfi mass and the unit 



m 



for 



charge. The calculations of the energy levels in this approximation are presented in 
different isotopes X of the hydrogen-like molecule {dtfi)Xee. Clearly, the treatment of dtfj, 
as a point-like particle is equivalent to the replacement of the exact potential V{\p + ar\) 
in the Schrodinger equation by the potential V{p) which describes the BO interaction 
between X and the point particle located at the dtfi center of mass. Thus, the effect of the 
dtfi structure, which leads to the shift of the zeroth-order energy levels, originates from the 
perturbation potential 

V, = Vi\p + ar\)-Vip) . (2) 

In the zeroth-order approximation Vp = 0, the solutions of Eq. with the total angular 
momentum L and its projection M are written as a product of the bispherical harmonics 
yt\ \Pi describing the angular dependence, the radial function of p describing the motion 
of nuclei in {dtp)Xee with the angular momentum /, and the radial function of r describing 
the internal motion in a mesic molecule with the angular momentum A. The unperturbed 
energies Eni and the corresponding square integrable radial functions $nz(p) of the {dtp)Xee 
vibrational and rotational states satisfy the equation 
f 1 



2p: 



2 



.1 9 + 



+ V-(p) - E,a 4.„,(p) = (3) 



dp \ dp J p^ 

where n is the vibrational quantum number. For the problem under consideration, one 
should consider both the bound and continuum states of the t + dp subsystem whose energies 
and wave functions satisfy the equation 

r 1 



2p 



Here E = —e^\ and 0(r) = (p^xi^) for the bound states and E = k'^/2pi and (f){r) = (pkxi^) 
for the continuum states with the wave number k. The functions (pvxi'f') are square integrable 
and the functions (pkxi'f') are normalized by the condition 

r^dr(j)l^{r)(j)gx{r) = 6{k - q) . (5) 







In correspondence with the Vesman mechanism, {dt^)Xee contains a dt^ mesic molecule 
in the weakly bound state with the binding energy Su {v — 1,\ — 1). Other dt^ states, 
whose binding energies significantly exceed all the characteristic energies of the problem 
under consideration, will not be taken into account in the calculation of the energy shifts. 

B. Perturbation theory 

The effect of the dt^ structure is small due to smallness of dt^ mesic molecule in com- 
parison with a characteristic length of dt^ motion in the molecular potential V{p). In other 
words, the perturbation Vp is small in comparison with V{p) and can be expanded in powers 
of the small parameter ar. Correspondingly, the dimensionless parameter of the perturba- 
tion theory is the ratio of the average distance between the dcuteron and the dt^ center of 
mass a{r) to the average amplitude of vibrations {p — a) in the molecular potential near the 
equilibrium internuclear distance a. 

One should note that the lowest-order term of the expansion Vp, which is proportional 
to ar, does not contribute to the energy shifts in the first-order PT; therefore, the energy 
shift of order (ar)^ must be obtained up to the second-order PT. Besides, Vp couples the 
rotational states with I = L ± 1 while the state with / = L remains uncoupled. As the 
separation of the rotational levels is comparatively small, the level coupling cannot be a 
priori neglected and requires explicit treatment. Thus, the energy shifts will be determined 
in the second-order degenerate PT by solving a secular equation 

detiV + + Sn- E]^0 (6) 

where and are the matrices with the matrix elements of the first- and second-order 
PT Vif^ and WH^, respectively, the matrix elements of Sn are {Eni + Su)Sii-^, and E is the 
level energy. 

The first-order PT matrix elements are 

V^ = / d'rd'pVp I Mr) r ^ni{p)^nH {P)yn''\p, r)yt'^{p, r) (7) 

and the second-order PT matrix elements include a sum and an integral over intermediate 
states describing simultaneous excitations of the exotic molecule with the quantum numbers 
V and I and a dtp mesic molecule with the continuum-state wave number k and the angular 
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momentum A 



wn 



E 



where 



(8) 



(9) 



Bearing in mind that the second-order PT calculation are of order a {r) , we expand Vp 
up to the second order in ar, which corresponds to the multipole expansion 

1 



dV 
dp 



2„2 



Vp = ar——Pi{cos6) H — a r 
f 6 



d'^V 2dV fd'^V ldV\ „ , 

+ -— + 2 I ^ - -— P2 cos( 



(10) 



dp'^ p dp \ dp^ p dp 

where 9 is the angle between the vectors r and p. The monopole and quadrupole terms in 
(I10|) contribute only in the first-order PT while the dipole one in the second-order PT. 

Calculation of the matrix elements VJp^ and W{1^ with the perturbation Vp in the form 
(fTUI) results in 

(11) 
(12) 



-a 



where 



' E UZU^n E A^ilUX)A^ihU\)hiE,e -E + Sn) 
[ux{k)]^ dk 



k^/2fii + A 



Q = J drr^ I ipii{r) ^ , 
ux{k) = / drr^(j)kx{r)(pn{r) 



M 

nl,nil 

Q 



p^dp^nl{p)^n^h{p) 



U^l,n,H='^Jp'dp<^nl{p)K,M 



D 

nLnili 



p^dp^M^mhip) 



'd^V 2dV\ 
dp^ p dp j 

fd^V _ ldV\ 
\ dp^ P dp j 

dV 



dp 



(13) 

(14) 
(15) 

(16) 
(17) 
(18) 



and the angular integrals A^{l\liXi) are given in the Appendix. 
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III. RESULTS OF CALCULATION 



A. Matrix elements 

The simple, though providing the required accuracy expressions for the multipole ma- 
trix elements (fTB|) . (fT7j) . and (fTH|) are obtained using the following reliable approximations. 



Firstly, the matrix elements are completely determined by the BO potential 



(U ppten 



for the hydrogen 



15|. As {dtfj.)Xee is 



molecule V{p) which is fairly well known from the calculations 
produced in low-energy t/i + DX collisions, only the lowest vibrational states should be taken 
into account. For these states are localized near the minimum of V{p) at the equilibrium 
internuclear distance a ~ 1.4 a. u., it is natural to use the harmonic approximation 

VH{p) = ^P2Uj\p-af + Vo (19) 

where only the frequency of vibrations uj is of importance for the calculation. Besides, an 
accuracy of the approximation ()19|) is estimated using the unharmonic approximation of the 
BO potential 

K(p) = ^Ai2tu2(p - a)2[l - «A/(p - a)] + Vo . (20) 

which takes into account the next term of the expansion in p — a. The approximation (j2l 
accurately reproduces the exact energies of the lowest vibrational states calculated in 

Secondly, the rotational energy in ^ for the hydrogen-like molecule Z(/ + l)/2/i2p^ ~ + 
l)/2/i2a^ ~ 10^^ is two orders of magnitude smaller than the vibrational energy u ~ 10^^. 
Therefore, under a usual approximation, the centrifugal term is treated perturbatively, i. e., 
the eigenenergies are given by 

Enl = EnO + Vrl{l + 1) (21) 

and the wave functions ^ni{p) will be taken independent of / in the same approximation. 
Indeed, the rotational spectrum calculated in jy, [D] is in good agreement with the above 
expression (PT|) with Vr ~ l/2/i2a^ ~ 10""^. Thus, under the above approximations, the 
radial wave function $ri«(p) i'^ potential p9|) coincides with the harmonic-oscillator wave 
function and the multipole matrix elements (HH), (fT7|) . (fTHj) are reduced to /-independent 
expressions 



UZ = \ ^ (v^5„_i,, + V^rTl5r.+,,^ , U::i = P2^', U^n = 2P2^' • (22) 



The unharmonic term of the potential Vu{p) leads only to modification of the dipole matrix 
element 



n+l,u 



- V 



(23) 



where the unharmonic correction is proportional to the dimensionless parameter t] = 
|aM/V2/i2tu ^ 0.14. 

Calculation of the quasi-nucleus matrix elements ()14l ITSj) is based on the smallness of 
the tfi size in comparison with the size of the loosely bound dtfi state {v = 1, A = 1). Thus, 
almost in all the configuration space tfi and d move as free particles and the bound-state 
wave function is approximated by 

^^11(0=^^6- (24) 

where k = y/2jlieu- The asymptotic expression has been widely used in description of 
the loosely bound states of mesic molecules and the asymptotic normalization constant Ca 
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171 . For the 



was determined by a comparison with the exact three-body calculations j9 , 
same reasons, the asymptotic expressions are used for the continuum wave functions (f)kx{f), 
viz., 

, , , [2 sinikr + 6o{k)) 
V TT r 



-k [cosS2{k)j2{kr) + sin S2ik)y2{kr)] 



(26) 



where 6x{k) are the tfi + d scattering phase shifts and j2{kr) and y2{kr) are the spherical 
Bessel functions. The d-wave phase shift S2{k) is actually very small, which allows either 

COS Jvl^ 

replacing y2{kr) by the leading term — of its asymptotic expansion or simply putting 

S^{k) = in Eq. Using the wave functions (j211), (123), and ^ with y2{kr) 

one obtains the quadrupole momentum 



cos kr 
kr 



Q 



5 CI 



8 



(27) 



and the expression 



4C2 
7r/ii£ii 



(28) 
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via the dimensionless integrals 



7 [sin^o(fc) - + 3)(A:/2) cos6o{k)f dk 
^'^'^ - J + 1)4(P + z) ' (29) 

7 [sin 62ik) + cos S2{k)f dk 
B. Shift and splitting of energy levels 

Energy shifts are obtained by solving the secular equation (jH)) which is reduced, due to 
the selection rules for angular momenta / and /i, to a 2 x 2 matrix equation for = L ± 1 
{L 7^ 0) and a scalar equation for I = li = L ^ 0. The energy shifts with respect to the 
unperturbed rovibrational energies Eni + eu are denoted as Ao(n) and A±{nl) for / = L 
and / = L ± 1, respectively. Note that the state with L = and / = /i = 1 is uncoupled; 
however, its energy shift A_|_(nl) will be determined in the same manner as for the other 
L 7^ states. 

The first-order PT matrix elements VJ"^ in Eq. © are calculated by substituting the 
radial integrals f/*^'*^ (j221), the quadrupole momentum Q (jTTj) . and the angular integrals 
A2{lllil) ()A4|) in Eq. (fTT|) . Note that V^"^ appears to be independent of the vibrational 
quantum number n and this index will be omitted in what follows. The matrix elements are 
scaled by a single dimensional parameter 

Vq = 7 r (31) 

16m2(mi + 7712 + m3)eu 
which is a characteristic energy for the problem. It should be mentioned that although 
Eq. ()3H1 does not contain a specific small parameter, vq turns out to be sufficiently small 
{vo/uj ~ 0.006) thus making the energy shifts small. As a result, one finds 

Vii^ = vo < (32) 
[25a,- Bu„ IJi = Lt1 

where the matrix elements Bu-, form the matrix 



1 2JL(L + 1) , 

B = — ^ , V ^ ) \ (33) 



in which the first row and column correspond to /, /i = L — 1 and the second ones to 
l,h = L + 1. 
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The second-order PT matrix elements Wi^_^ in Eq. (jH)) are calculated by substituting 
Uji (US) and /a (A) ^ in Eq. which gives the expression 

Wll^ = -vo^ E Y.^axin,u)A^{mX)A^ihUX)Jx + (34) 

via the energy scale vq and the dimensionless factors. Solving the secular equation ©, 
one can safely replace, up to an accuracy of the second-order PT, the eigenvalue E in the 
argument of Jx by the unperturbed value Eni. Thus, the calculation of the energy shifts is 
basically accomplished by derivation of Eqs. ()32ll34j) . 

However, it is reasonable to make further simplification of ()34|) by neglecting the difference 
of the rotational energies in the argument of Jx, which allows obtaining an explicit and 
sufficiently accurate dependence of the energy shifts on the quantum numbers n and /. As 
the rotational energy is much smaller than the vibrational quantum uj, one replaces the 
energy differences En±ii — E in the argument of Jx by the /-independent values -Eniio ~ 
Eno = ±u!. Using the angular integrals Af{llli\) ()A5l IA6|) and introducing the notation 
J^ = Ja (1 ± tu/eii) for integrals independent of n and / one obtains 

Wil=Vo{ (35) 
2 - Pn)5u, + a„ - l)Bu, , l,h = LTl 



where 



IQuj 

a„ = 1 - 



[n + 1){J+ + ij2+) + n{J^ + i J2" 
[n + l)(Jo+ + ht) + n{J, + 



(36) 
(37) 



Siren 

determine the explicit dependence on the vibrational quantum number n. As a result, the 
sum of Vii^ (jH^ and Wll_^ takes a simple form 

T/ MU,n l = h = L (L ^ 0) 

Vih + =vo< , (38) 

[ Pn^ih - anBii^, /, /i = L =F 1 

i. e., the parameter /3„ determines the constant shift VoPn of all level energies Eni whereas a„ 
determines the level splitting. Using and (j7T|) in the secular equation ^ one obtains 

Ao(n) = VoiPn - an) , (39) 



A±(n/) = VoPn T VMI T 1) + 1] ± V^O«n + '^VoanVr + [2(/ T I) + l?V^r ■ (40) 
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The effect of coupling of the rotational states with / = L ± 1 is explicitly taken into account 
in expression (^UJ. Generally, the effect decreases with decreasing ratio of the level splitting 



to the energy difference between the rotational states 



voan 



i. e., with increasing total 



Vr{2L + l) 

angular momentum L. As follows from the numerical values of Vq, Vr, and a„ for all the 
considered states n = 2,3 (Section IIII (]|) . even in the worst case L = 1 the energy shifts 
calculated with and without allowance for the coupling of the rotational states differ at the 
most by O.OlmeV for {dtfi)dee and O.OSmeV for {dtfi)tee. As these values are beyond the 
accuracy of the present calculation, it is quite reasonable to neglect coupling, i. e., to use 
the diagonal approximation for the secular equation (jH}, which allows obtaining a simple 
expression 



A±(n/) = Vq 



1 



^^0 



Ctr. 



2L 



(41) 



2{lTl) 

Note that Eqs. (HOI) and (jUI) are valid both for L = and / = 1 when A+{nl) = 
VoiPn + «n) and for L = 1 and / = when A_{nO) = fo(/?„ — ^n/S). The sign of a„ 
determines the relative position of the levels so that the energies satisfy the inequalities 
EnL+i > EnL-1 > EnL foT > and the inverse inequalities for a„ < 0. The largest 
energy splitting is predicted for / = 1 between the states with L = and L = 1, viz., 
A+(?2l) - Ao(n) = 2voan. 



C. Numerical results 



The energy shifts and level splittings will be calculated by solving the eigenvalue equa- 
tion © using formulas (jHU-EH)- In addition to the particle masses = 206.768a.u., 
rrid = 3670. 484a. u., and rrit = 5497.922a.u., calculation of the matrix elements VJ/^ and W^^ 
requires the vibrational u and rotational Vr energies of the exotic molecule {dtfi)Xee, the 
binding energy eu and the asymptotic constant Ca of the dtfi loosely bound state, and the 
low-energy tfi + d scattering phase shifts Sx{k) which determine the integrals Jx (l + "^^7^)- 

The vibrational quantum cu and the rotational-energy constant Vr are determined by the 
BO internuclear potential of the hydrogen molecule near its minimum or, equivalently, by 
the low-lying part of the {dtfi)Xee vibrational-rotational spectra calculated in j^I^. Fitting 
the BO potential near the equilibrium distance a = 1.401 to the harmonic, unharmonic, and 
Morse potentials provides consistent determination of both lv and the parameter aM- As the 
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BO potential is independent of the isotopic composition, both /i2Co'^ and are independent 
of the masses of heavy particles due to Eq. (j^UI) . The result of the fit gives uj = 321.8meV 
for {dtfi)dee (correspondingly, u = 273.1meV for {dtfi)tee) with a few per cent accuracy and 
the parameter om = 0.7. For these parameters, the energies of the lowest vibrational states 
in the approximate potential are in reasonable agreement with the results of lal?!. 

fin 

The rotational spectra calculated in ja l2| are fitted to Eq. for 1 < / < 10 and each 
1 < n < 4. For the lowest vibrational state n = 1, one obtains Vr = 2.43meV for {dtiJ,)dee 
and Vr = 1.85meV for {dtfi)tee. These values agree with the simple estimate 2^2Vr ~ 1/a^ 
that determines the isotopic dependence of Vr- Although Vr slightly decreases for the higher 
vibrational states, the above values will be used for n > 1, which leads to a few per cent 
error. 

Determination of the binding energy eu of a dtfi loosely bound state was a subject of 
numerous elaborated calculations. As a result, the value en = 596meV [l|, 0,Q is obtained 
for the lowest hyperfine state by taking into account relativistic effects, hyperfine effects, 
finite nuclear size, and vacuum polarization. The asymptotic constant Ca was determined 
in a number of papers P,^, 12 1 by a comparison of the asymptotic expression ()24|1 with the 



three-body wave function. In the following, it is accepted the value Ca = 0.874/v^ obtained 
in the latest elaborated calculation ^ of the wave function in a wide asymptotic region of 
large distances between d and tfi. Using Ca, eu, and u one can calculate the energy scale 

vq (jSD). As cu^ ~ l//i2 = — 1^ the parameter vq dSU is independent of m^, i. e., it 

(mi + m2)m3 

is the same for any isotope X. Given the above numerical values one obtains Vq = l.SlmeV. 

For the sake of completeness, it is interesting to estimate the energy scale for the molecule 
{ddjjL)Xee too by using the values en = 1975meV, Ca = 1.006/-\/2, and uj = 257meV, which 
gives Vq = 0.4meV. Although the present approach requires some modifications to describe 
{ddn)Xee, viz., taking into account the identity of nuclei in ddfi and the essential role of the 
unharmonic corrections to the BO potential, one can qualitatively conclude that the energy 
shifts in {dd^)Xee are 4 — 5 times smaller than in {dtfi)Xee. 

In the present approach the energy shifts in the first-order PT are given by simple de- 
pendence on the angular momentum / flH^ - IH^ containing a single parameter vq. It is 
worthwhile to compare this result with the first ever elaborate six-body calculation of the 
{dtiJ,)dee energy shifts in the first-order PT 8]. In this paper, the molecular structure, i. e., 
the dependence on /, was explicitly taken into account in contrast with previous calcula- 
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tions |9|,lld,ll8| where the /-independent energy shift was obtained by scahng the result for 
the atom-hke four-body system {dtfi)e. As pointed out in this paper, the monopole con- 
tribution calculated in Ref. O] depends on the choice of the coordinate system that does 
not allow a comparison. For this reason, only the quadrupole contribution to the first-order 
PT energy shifts of Ref. 0] will be compared with the present results. The quadrupole 
contribution of Ref. ^ is given in Table II of that paper, while in the present approach 
AE^Q^ = {10/3)voA^{llll), as follows from Eqs. (HD), and (EZj). The both resuhs are in 
excellent agreement with each other, as shown in Table H] Note that in the present approach 
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2 -1/5 


-1.21 


-1.17 


1.76 
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2 4/25 


0.96 


0.94 


1.77 



TABLE I: Quadrupole contributions Ai^^ (in meV) to the first-order PT energy shifts of the 
present calculation and those from Ref. [g] for different / and L. Also presented are the angular 
integrals A2{llll) and the parameter vq corresponding to the energy shifts of Ref. ^. 

the dependence on angular momenta is completely determined by the factor ^2 (/l/l) which 
is also presented in Table HI To a good accuracy, the results of Ref. reveal the same 
dependence on angular momenta which approves the description of energy shifts by a single 
parameter vq. To emphasize this fact, the quadrupole correction calculated in Ref. ^ is 
expressed in the form AEq^ = {10/3)voA2{llll) with the variable Vq presented in Table El 
Indeed, vq is practically independent of / and L and agrees with vq = l.SlmeV. Agreement 
between the present one-parameter result for the quadrupole correction and the elaborate 

n 

six-body calculation |8] is a good argument for the validity of the present approach. 

At last, one should obtain W^^ which requires evaluation of Jx (l + ^^f^^ by 

using the tji + d scattering phase shifts 5x{k) in the integrands of Eqs. (^^1) and (jHUI) . The 
■ow-energy sc attering phase shifts were determined in a number of three-body calculations 



IS 



2G 
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24| . whose results are in good agreement with each other. Using 5\{k) from 
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these calculations and integrating and (I^Uj) in the energy interval < k"^ /2jj,i < lOeV, 
one obtains J\ {l + ^^f^^ with a relative accuracy about 0.01. 

Calculating the matrix elements Va-^ (jH^ and W[l^ (jH^ and solving the eigenvalue equa- 
tion (jni) one obtains energy shifts presented in Table |n] for {dtfi)dee and {dtfi)tee. Note 
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1.09 
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TABLE II: Energy shifts (meV) for a few states of {dtfi)dee and {dtfi)tee with the vibrational 
quantum number n, the total angular momentum L, and angular momentum / of the hydrogen- 
like molecule with the point-like dtfi quasi-nucleus. 

that applicability of the harmonic approximation for the BO potential was checked by using 
the modified dipole matrix element (f^Hj) in the calculation, which gives an estimate of the 
unharmonic correction of the order of 5% in the energy shifts. Calculations reveal that 
the energy shifts are essentially dependent on the isotopic composition and the molecular 
quantum numbers n and /, which is basically connected with the cancellation of the first- 
and second-order PT contributions. In particular, the energy shifts decrease with increasing 
n so that A± become very small or even negative for n = 4. The reason for this dependence 
is an increasing in the dipole matrix element (|^^ with increasing n, which, in turn, leads 
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to an increasinEr in the second-order PT contribution. The cancellation effect was widely 
discu..d. e. inflHEi; „eve.Me.. t.e dependence on .he n.o,ecu,a. .nantun. .ate 

was beyond the scope of those papers^here only the atom-like system {dtfi)e was calculated. 
On the other hand, the calculation 1,8.] determined the /-dependence only in the first-order 
PT. 

The dependence of the energy shifts on quantum numbers is illustrated in Figure ^ for 
the {dtii)dee states with n = 2, 3 and / = — 4. In addition to a decreasing in the energy 

2 I , , ^ 

A ^ 



1.2 - 



0.8 



0.4 I ' ' ' J 

12 3 4 

1 

FIG. 1: Energy shifts (meV) of {dt^)dee. Crosses, squares, and circles denote, respectively, A+, 
A_, and Aq. The results obtained by simplified formulas H41() are denoted by the solid, 

dashed, and dotted lines, respectively. Three upper curves correspond to the vibrational quantum 
number n = 2 and the lower ones to n = 3. 

shifts for higher n, notice the inverse ordering of levels, i. e., the highest level with L = 1 — 1 
for n = 2 becomes the lowest for n = 3. Except for A_|_, whose values at small / = 1, 2 are 
quite different, the results reveal weak dependence on I with splitting of levels of the order 
of 0.2meV. 

As discussed at the end of Section lTlI Bl the dependence of energy shifts and level splitting 
on quantum numbers are expressed to a good accuracy by simple formulas (jSHl), (ED) via 
few parameters. Numerical values of and J^2 given in Table UTTl for {dtfj.)dee and 
{dtii)tee. As is clearly seen in FigureHJ the simplified expressions (jHS)), provide a reliable 
description of energy shifts. 
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d 0.917 0.113 0.236 0.051 0.073 
t 0.778 0.117 0.216 0.052 0.069 

TABLE III: Dimensionless parameters entering into the simplified expressions (|39)1 . (|41|) for the 
energy shifts of {dtfi)Xee. 

IV. DISCUSSION AND CONCLUSIONS 

The hydrogen-like molecule {dtfi)Xee is treated within the framework of the three-body 
model for heavy particles d, tfi, and X. The model is based on the fact that a tfi mesic atom 
is small in comparison with its mean separation from a deuteron and that the size of a dtn 
mesic molecule is small in comparison with the amplitude of vibrations in {dtfi)Xee. In this 
approach, the interaction of the charged particles d and X is described by the well-known 
BO potential of the hydrogen molecule, while the description of a dtfi mesic molecule is 
given in terms of the binding energy ^n, the asymptotic constant Ca, and the low-energy s- 
and (i-wave scattering phase shifts 6x{k) regardless of the explicit form of the d + tfi effective 
potential. 

In the present approach, the shift and splitting of the {dt^)Xee energy levels which 
result from the internal structure and motion of the dtfi quasi-nucleus are calculated in 
the second-order PT. This allows one to find the energy levels, i. e., the positions of the 
tfi + DX scattering resonances with an accuracy about a tenth of a meV, which is of key 
importance for determination of the dtfi formation rate. Calculations are performed for 
different vibrational (n = 2, 3) and rotational (0 < / < 4) states for the molecules of the 
different isotope composition X = d,t. In this respect, note that different vibrational states 



23- 



It should be 



of {dtjj,)Xee can be currently observed in the atomic beam experiments 
emphasized that the effect of the dtfi structure removes the degeneracy of unperturbed states 
with the same I and different L, which produces a triple-resonance structure in place of every 
unperturbed level except the one with / = 0. As the splitting value is of the order of the 
shift itself, the effect of splitting should be taken into account in the energy dependence of 
the resonance formation rate. 

The following aspects of the present calculation are worth mentioning. The first-order 
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PT quadrupole contribution to the energy shifts is in agreement with the elaborate six- 



body calculation |8[, which is a good argument for the validity of the present approach. 
Furthermore, it is shown that for all the considered states the effect of coupling of the 
rotational states with / = L — 1 and Z = L + 1 is beyond the accuracy of the present 
calculation. In addition, the energy shift and splitting is conveniently expressed by simple 
analytical dependence (jSH), ()41|) on the rotational- vibrational quantum numbers n and /. 

It is of interest to compare the present result with the direct non-perturbative three-body 
calculation of tjji + D2 scattering (lll.ll2| in which the effect of the dtfj, structure is explicitly 
taken into account. The positions of narrow resonances calculated in this paper correspond 
to energy levels of {dt^)dee for L = (/ = 1) and n = 3,4. For two approximations of the 



effective potential between tfi and d used in 



11 



12| . energy shifts are, respectively, ImeV 



and 4meV for n = 3 and 2meV and 5meV for n = 4. The cause of the noticeable difference 
(about 3meV) is not clear since both potentials allow a good description of the low-energy 
properties of the + d system. The dependence on the choice of the effective potential 
and a limitation only by L = hinders a quantitative comparison of the present results and 

11 . 1^ exceed the present 



those of [ill, • Qualitatively, the energy shifts obtained in Ref. 
ones and, contrary to PT considerations, the value for n = 4 is higher than for n = 3. This 
n-dependence clearly deserves further investigation. 

Finally, it should be mentioned that the present approach, which reliably takes into ac- 
count the structure of the exotic molecule, is promising for wider applications, in particular, 
for determination of resonance positions and formation rates beyond PT by solving the scat- 



3, 



tering problem. Till now, except Refs. [UJ, ll3 , the formation rates have been calculated 
only in the first-order PT. In this respect, the result of Ref. j26| shows that the first-order 
PT dipole approximation is questionable and one should do more refined calculations. In 
addition, it is of interest to apply the present a ppr oach to the problem of the resonance 



formation of metastable dt/i mesic molecules 
with D2 molecules. 







l28l | in collisions of excited t/i mesic atoms 



APPENDIX A: ANGULAR INTEGRALS 

The following angular integrals are necessary to calculate the matrix elements 

AUl^i) = J dpdiPK{cosd)yl^^''\p, r)y[;ll{p, r) (Al) 
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where 9 is the angle between two unit vectors f = r/r and p = p/ Pk{x) is the Legendre 
polynomial, and the bispherical harmonics are defined as 



(A2) 



Evaluating the integral ()A1|) one comes to the expression in terms of the Clebsh-Gordon 
coefficients and 6j-symbols 



\h+L 



{21 + 1)(2A + 1){10K0 I liO){XOKO \ AiO) 



k I K 
X Xi L 



(A3) 



The matrix elements Vif^ flll|) are expressed in terms of the integrals ()A3|) with K = 2 
and A = Ai = 1 which are explicitly written as 



3 JL{L + 1) 



(A4) 



5 (2L + 1) 

The matrix elements W[l^ (|T2|l are expressed in terms of the integrals ()A3|1 with K = 1, 
A = 1, and either Ai = and /i = L or Ai = 2 and I = L ±1. The explicit expressions read 

1 



Af(/1L0) 



Af(/l/i2) 



3(2L + r 



1 

I 2\ 



(L + / + 3)(L + / + 4)(L - / - 2)(L - / - 3) 
15(2Z + l)(2/ + 3) ^ 



(L - / + 2)(L - / + 3)(L + / - 2)(L + / - 1^ 
15(4/2 - 1) 



h = l + l 



h = l-l 



(A5) 



(A6) 
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